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Abstract 

Experiments determining the lifetime of excited electrons in crystalline copper 

reveal states which cannot be interpreted as Bloch states [S. Ogawa et al, 

Phys. Rev. B 55, 10869 (1997)]. In this article we propose a model which 

explains these states as transient excitonic states in metals. The physical 

background of transient excitons is the finite time a system needs to react to an 

external perturbation, in other words, the time which is needed to build up a 

polarization cloud. This process can be probed with modern ultra-short laser 

pulses. We calculate the time-dependent density-response function within 

the jellium model and for real Cu. From this knowledge it is possible within 

linear response theory to calculate the time needed to screen a positive charge 

distribution and - on top of this - to determine excitonic binding energies. 

Our results lead to the interpretation of the experimentally detected states as 

transient excitonic states. 

78.47.+p, 72.15.Lh, 71.10.-w, 71.20.Gj 
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I. INTRODUCTION 

Measurements of the lifetime of excited electrons in CuEB using time- resolved two-photon 
photoemission (TR-2PPE) spectroscopy! in connection with very short laser pulses reveal a 
couple of amazing results. Most stunning is the fact that the data show lifetimes of states 
along the (100) and (111) direction (T — X and T — L, respectively) with energies between 
1 and 4 eV above the Fermi energy Ep. As can be seen from the band structure shown in 
Fig. [I], Cu has no Bloch states at these energies in these two directions. But also the data for 
the lifetime of electrons in the remaining direction (110) - the T — K direction - are puzzling. 
Ab initio calculations of the lifetime of hot electrons along the (110) direction using the GW 
approximation!'! reveal a good agreement with the experimental data for electron energies 
between 2 and 3 eV. For energies lower than 2 eV the results of the GW calculation deviate 
markedly from the experimental data. This is particularly remarkable because typically 
many-body calculations within the GW approximation agree very well with experimental 
data. This holds for both, the real part of the resulting quasiparticle energies!'! as well as its 
imaginary part,!&!il respectively. The imaginary part of the quasiparticle energy of a band 
state is inverse proportional to the lifetime of this state.0 Therefore the discrepancy of the 
theoretical and experimental results must be due to physics which is not covered by the GW 
approximation. The last point about the experimental data we would like to mention is the 
fact that the lifetimes in the three directions are very similar! although the band structure 
gives no hint for this experimental finding. 

The discussion in the preceding paragraph shows that TR-2PPE experiments detect 
states which are obviously no band states (which are single-particle states). Consequently 
there has been some speculation in the literature about the origin of these states, most of 
which violate well established fundamental physical laws. The only explanation which is 
consistent in this respect and does involve only single-particle states is that electrons which 
originally had a momentum in the, e.g., (110) direction are scattered from a (110) low-energy 
electron diffraction (LEED) state into a, say, (100) LEED state. This would explain why 
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experiments measure the lifetime in this direction. However, the cross section for this process 
should be small. Moreover, this mechanism would not explain the difference between the 
results of calculated lifetimes within the GW approximation and the experimental ones. In 
order to escape these difficulties Cao et al. speculated^ that "the strongly localized 3d holes 
generated by photoexcitation of 3d electrons can trap excited electrons through attractive 
Coulomb interaction". We picked up this suggestion and outlined the physical mechanism 
of excitons in metals on a very short time scale, so-called transient excitons.0 In this article 
we present a model of how to actually calculate these states. 

Now it is common wisdom that excitons exist in semiconductors and rare gases but not 
in metals.ll! However, in these cases one is thinking of stationary conditions under which the 
valence electrons of the system have enough time to built up a polarization cloud. For large 
times the initial Coulomb interaction between the hole and the electron transforms into a 
Thomas-Fermi-like potential in metals0 and a Coulomb potential which is screened by the 
macroscopic dielectric constant in semiconductor^. But this polarization needs a final time 
to be built up and it is during this time that transient excitons can exist .0 Consequently 
transient excitons cannot be detected with traditional photoabsorption spectroscopy. It 
needs modern TR-2PPE spectroscopy combined with the ultra-short laser pulses which are 
available nowadays. 

The article is organized as follows. In Sec. |H| we present a model that allows us to 
calculate - and therefore proof the existence of - transient excitonic states. In the two 
following sections we apply this model to the jellium model and to real copper. We finish 
with some estimates about the lifetime of electrons in transient excitonic states. 

II. A MODEL FOR TRANSIENT EXCITONS IN METALS 

In this section we present a model for the calculation of transient excitonic states. A com- 
plete treatment of the subject implies the solution of the Bethe-Salpeter equation 
using the time-dependent particle-hole potential v e h(t) as interaction. However, it is not our 
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objective to perform an ah initio calculation in order to calculate the excitonic energies to 
high accuracy. The intent of this article is to show that transient excitonic states in metals 
do exist. For this, a simple but nevertheless realistic model is very useful. 

Starting point of our discussion is the potential due to a localized hole in a 3d band. We 
obtain it from the density of a 3d-band electron as 

Ja|^(r)| 2 r G MT 
n h {v) =1 , (1) 

1 else 

where the Bloch function fzdir) is the solution of the Kohn-Sham (KS) equations^ 

h 2 

[ - T^V 2 + Veff(r)] </? qj -(r) = e qj </? qj (r) . (2) 

q and j denote a wave vector in the Brillouin zone (BZ) and a band index, respectively. 
v e s(r) is the mean-field potential in which the KS electrons move. Its exchange-correlation 
potential is treated within the local density approximation (LDA).0 The wave functions 
Vqj( r ) are expanded with respect to plane waves and the interaction between the ionic cores 
and the valence electrons is described by an ah initio norm-conserving pseudopotential.0iil 
a has to be chosen in such a way that 



d 3 rn h (r). (3) 



MT 



In Cu we use for the radius of the sphere MT a value of 2.3 Bohr. The 3o?-hole density n h (r) 
is expanded with respect to spherical harmonics and using Poisson's equation the potential 
corresponding to n h (r) is obtained as 

l,m ^ ^> 

nf m (r) are the (l,m) components of n h (r). r< and r > denote the smaller and larger length 
of r and r', respectively, f represents the angular part of r. Figure shows the first non-zero 
components of v h (r). The by far dominating component is the spherical symmetric part. 
In the following we approximate v h (r) by its spherical part and treat this as an external 
potential switched on at t = 0, 
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v cxt (r, t) = 4vre / dr> r' 2 — n %(r') K , (r) Q(t) , (5) 

perturbing the system. 

The response of the system to the perturbation v ex t(r, t) is calculated within linear re- 
sponse theory (LRT).H In LRT the induced density due to the perturbing potential is given 
by 

rw(r, t) = / dV / di! X (r, r'; t - f) v ext (r', t') , (6) 

where x( r , r '; t — t') is the (retarded) density-response function. Since we are dealing with a 
system with translation symmetry it is useful to transform Eq. (|6|) to Fourier space, 

n- md (q+ G,uj) = J2 XG,G'(q,^) VertCq + G'jW) . (7) 



G' 



XG,G'(q5 ^0 i s obtained from the polarizabilityi! 

< Me- i(q+G) *|k + q,j' > < k + q,jV (q+G ' )x |k,j > (8) 
by solving the matrix equation^! 



XG,G'(q> w ) = PG,G>{<l,Uj) 



J2 Pg,Gx(cI,^) 

Gi,G2 



S Gl , G2 « Gl (q) + r c (G!-G 2 ) 



XG a ,G'(q> w ) ( 9 ) 



for each wave vector q and frequency cu, respectively. In Eq. (§) the f^j denote the occupa- 
tion numbers, e^j the eigenvalues, and < r|k, j >= <^k,j(r) the wave functions as obtained 
from solving the KS equations ([J). V is the normalization volume. The sums run over all 
wave vectors k in the BZ and the band indices j and j', respectively. In Eq. (|9|) Ug(q) is 
the Coulomb potential, t>G(q) = 47re 2 /|q+ G| 2 , and / XC (G) represents a vertex correction 
connected to the local-field factor G by / XC (G) = -v G (q) G G (q). Within the LDA / XC (G) 
depends only on the reciprocal lattice vector and not on the wave vector. Setting / XC (G) = 
results in the random phase approximation (RPA) for \. The time-dependent local density 



approximation (TDLDA) is given by / XC (G) = / d 3 r e~ tGr dv^ c (r)/dn(r), where v xc (r) is the 
exchange-correlation potential used in the KS equations (|2|). 

Having calculated the induced density n in d(q+G, lj), the induced potential Vi n d(q+G, uj) 
is obtained via Poisson's equation. Adding the external potential f ex t(q+ G,uj) and using 
the Kramers-Kronig relation,0 a Fourier transform with respect to frequency yields the 
following expression for the total potential, 

Utot(q + G,t) = Y^ ^ext(q + G') x 

G' 

Attf 2 1 r°° 1 / \ 1 

8 GtG > + i — ^ - / du-lm XG M<l>u) H 6(t) , (10) 

|q + Lw\ z 7T J-oo C<J V / J 

where t> ex t(q + G) is the frequency-independent part of f ex t(q + G,uj). A final Fourier 
transform leads to the total potential as a function of r and t, 

i BZ 

«w(r,*) = ^EE ei(q+G)r ^ot(q+G,t) . (11) 
q G 



III. THE JELLIUM MODEL 

Before we treat real Cu with our model, let us first apply it to the jellium model. Even 
here it leads to a couple of interesting results. In the jellium model, where there are no d 
bands, we use for the external potential a bare Coulomb potential 

v est (r,t) = --Q(t), (12) 
r 

i.e., we calculate the time-dependent screening of a suddenly created point charge in an 

electron gas.0 Because of the isotropy assumed in the jellium model the potential, polariz- 

ability, and density-response function are simple scalar quantities. Equation (|) reduces to 

the Lindhard functionB and Bq. © becomes a scalar equation. So Bq. © simplifies to 

47re 2 f Aire 2 2 r°° 1 / \ 1 

vtat(q,t) = —2- 1 + — r - / du-lm X {q,uj) 1 - cos(cut) 9(f), (13) 
q { q 7r Jo u V /J 

where we used the fact that lmx{q, oj) is an odd function with respect to uj. Consequently 

Eq. (PI) becomes 



Vtot(r, t) = J d 3 q e^ r v tot (q, t) . (14) 

In the jellium model Imx{q, oj) has a delta-function-like shape for small wave vectors q, 
representing the plasmon excitation. It is centered at the plasmon frequency u p with uo 2 = 
Aime 2 /m, where n is the electronic density. In order to handle this structure a numerical 
broadening rj has to be introduced. If not said otherwise, we use a broadening of rj = 0.01 
eV, a value corresponding to the real width of the plasmon excitation in alkalines.il In order 
to be sure to sample the plasmon peak properly we used a frequency mesh with a spacing 
of 0.5 meV (and less) for wave vectors up to 1.5 /cp. 

The upper panel of Fig. |3] shows a plot of v tot (q, t) as a function of t for a Wigner-Seitz 
radius of r s = 2.07 and a wave vector of q — 0.1 /cf, hp being the Fermi wave vector, x was 
calculated within the RPA. The figure displays Vt Q t{q, t) for three different values of rj in a 
time window of just one femtosecond. The solid, dotted, and dashed curve correspond to 
values of rj of 0.01 eV, 0.1 eV, and 0.5 eV, respectively. Even on this short time scale the 
influence of r\ can be noticed; v tot (q, t) decays at different rates. The frequency of v tot (q, t) 
for the three values is nevertheless the same. It is the plasmon frequency u p . For r s = 2.07 it 
is uj p = 15.83 eV which corresponds to a period of 0.26 fs. That the potential for small wave 
vectors oscillates with the plasmon frequency can also be seen if lmx{q, oj) is approximated 
by the plasmon-pole approximation (PPA) 

limlmx(g,^) = -^-u p [d(uj - u p ) - 5{u + u p )\ , (15) 

where v q = 47re 2 /q 2 . In this case the frequency integral in Eq. (|13|) can be done analytically 
leading to 

Attp 2 

vZ\q,t) = — r Q{t)^ p t). (16) 

In the case of small wave vectors the total potential should approach the Thomas-Fermi 
potential 



for large times, ^tf is the Thomas- Fermi wave vector, q^ F = Qirne 2 



j Ep. Note that Eq. flTTD 



is derived for the limiting case of q — > and uo —>■ 0. Therefore it cannot be expected that 
lim t ^oo ftot(<Z) t) = v TF (q) holds for large wave vectors. For small wave vectors it is fulfilled to 
a high degree. As an example consider v to t(q, t) for q = 0.1 kp and t — > oo. For r/ = 0.01 eV 
the difference to the Thomas-Fermi value of 10.57 a.u. is just 0.8%. The deviation increases 
for larger broadenings. 

Having discussed the total potential as a function of wave vector we now turn to ftot(^, t). 
In order to perform a good quadrature of the integral in Eq. (|14]) we calculated v to t{<l,t) on 
a wave vector mesh which is dense enough to sample the oscillations of Vtot{<l, t) with respect 
to q. They are shown in the lower panel of Fig. [| where we plot v tot (q, t) as a function of wave 
vector at an intermediate time. In general the frequency of the oscillations increases with 
increasing t and with increasing q. At the same time the amplitude of the oscillations die 
out quicker so that at large times the asymptotic potential is approached. The plot shows 
another feature which holds for all times; the oscillations stop between 0.7 and 0.8 kp, i.e., 
they stop in the wave vector range where the plasmon excitation is destroyed by particle- 
hole transitions (Landau damping) and the density-response function becomes dominated 
by single-particle transitions. 

Figure f| shows snapshots of v t ot{f,t) as a function of r for fixed times calculated for 
r s = 5 (upper row) and r s = 2.07 (lower row). At t — 0, v tot (r, t) is just the Coulomb 
potential (upper left plot). The horizontal lines are the first four eigenvalues for this potential 
calculated using a reduced mass of /i = 0.5 due to the equal masses of the hole and the 
electron, respectively. The next plot shows v tot (r, t) at a larger time. The potential is already 
narrower and consequently the eigenvalues are shifted upwards. This trend continues and 
at t = 0.25 fs v tot (r,t) crosses the Thomas- Fermi potential [Fourier transform of Eq. (|Hj)1 



(dashed line) for the first time. At this time the potential has no bound states any more. 
The system is now overscreening the perturbation which leads to a decaying oscillation of 




r) = e 

r 



-<?TF r 



(18) 
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Vtot(r, t) around its asymptotic form which lasts for 10 to 20 fs, the exact value depending 
on r s . This is shown for r s = 2.07 in the lower panel of Fig. |j. At t — 20 fs the potential 
has reached its asymptotic value. As can be seen from the last plot, the difference between 
this potential and the Thomas-Fermi potential is rather small. 

Before we close the discussion about the jellium model we want to discuss the influence 
of the vertex correction / xc . We discuss this for r s = 5 since / xc for this Wigner-Seitz 
radius is more than eight times larger than the f xc for r s = 2.07 (-24.443 a.u. and -3.917 
a.u., respectively). We use the parametrization of the exchange-correlation potential due to 
Perdew and Zunger.il The upper panel of Fig. [5] shows a comparison of v to t(q,t) calculated 
with (dot-dashed line) and without (solid line) the use of vertex corrections, respectively. 
The figure shows that for wave vectors smaller than 2.2 kp the two curves start to differ. 
For small wave vectors the difference is quite remarkable. After the Fourier transform the 
difference in the resulting Vx t{r, t) is much smaller. This is because the difference for small 
q in v tot (q,t) is reflected in a difference for large r in t> tot (r, t). For our considerations this 
range is not of importance. 



IV. TRANSIENT EXCITONS IN COPPER 

The heart of the calculation of the total potential Eq. (10) is the calculation of the 
polarizability given by Eq. @ for all wave vectors in the BZ and a sufficient number of G 
vectors needed to obtain a reliable density- response function Eq. (|9]). In Fig. | we show the 
imaginary part of xg,g , ( ( 1- u; ) for G = G' and three wave vectors q + G. The upper panel 
of Fig. ^ shows the imaginary part of the density-response function for a small wave vector. 
It displays clearly the characteristic double peak structured around 20 and 29 eV sitting on 
top of a broad structureless hump which decays only very slowly. The calculations of the 
polarizability were done using 200 bands, i.e., considering transitions up to 355 - 390 eV 
above the Fermi energy (the exact number depends on the wave vector). This energy cutoff 
is the reason for the sudden decline of Imx G ,G'(q> u ) between 360 and 380 eV which can be 
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observed in all of the three plots. The middle panel shows Imxcd^i, w) for a wave vector on 
the edge of the BZ and the lower panel the imaginary part of the density-response function 
for a wave vector q+ G in the third zone. In all cases XG,G'(q> w) was calculated within the 
RPA using a numerical broadening of rj = 0.4 eV in the calculation of the polarizability. It 
is obvious from the plots that this value is small compared to the broad structures of the 
density-response functions of Cu. 

In the jellium model the screening process was completely determined by the plasmon 
excitation. The plots in Fig. ^ make clear that this is not the case for Cu. The density- 
response function of Cu is dominated by single-particle transitions and the already mentioned 
dispersionless double-peak structure caused by collective excitations of the d electrons.il 

As was already mentioned, the external potential Eq. ([|) is spherical symmetric. Al- 
though the full matrix of Xg,G'(cIj ^) wa s used in the calculation of the total potential 
^tot(q + G,t) (we used the first two shells of reciprocal lattice vectors, i.e., 15 G vectors) 
this leads to a nearly spherical symmetric total potential. Therefore the Fourier transform 
to real space could be done replacing the sum in Eq. ( |TTD by an integral and performing the 
angular part of the integration analytically. Figure [F] shows the total potential in Cu as a 
function of r for four different times. The plot in the upper left corner shows the potential 
at t = fs, i.e., the external potential. For finite times the potentials quickly becomes shal- 
lower. This behavior is the pendant to the time dependence within the jellium model where 
the potential became narrower. As in jellium model the total potential is overscreening 
(see lower left panel in Fig. [j]) and oscillating around the equilibrium potential. However, 
compared to the jellium model this process is faster in Cu. The oscillations are completely 
died out at 2 fs. At this time i>tot(?", t) has reached its asymptotic value. 

We are now in the position to apply the textbook approach for the estimate of excitonic 
binding energies in semiconductors to copper. We chose a high lying d band and (typically) 
the lowest unoccupied band for a given direction in the BZ and calculated the effective 
masses for each of the bands (m* h and m* for the effective mass of the hole in the d band and 
the effective mass of the electron in the conduction band, respectively). A typical case is 
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given by the two bands drawn as thick lines in the L — T panel of Fig. |8|. From the knowledge 
of m* h and m* the reduced mass \i for each wave vector is determined. Then the Schrodinger 
equation for t>tot(?", t) at a given time t is solved for each /i. t has to be chosen in such a way 
that v to t (r, t) still possesses bound states. The resulting lowest eigenvalue is plotted with 
respect to the eigenvalue of the unoccupied band. This is the excitonic binding energy. 

In Fig. [| the excitonic binding energies in the three main crystallographic directions are 
drawn as horizontal bars. The bands which were used in the calculation are indicated by 
thick solid lines. Since in the T — K direction there are two unoccupied bands close together, 
the calculation was done for both bands. Figure |8] shows that excitonic states in the L — T 
direction exist in an energy range where there are no band states. The energetic positions 
of the excitonic states are in agreement with the experimental data (see Ref. p]). The 
excitonic energy levels were determined from the eigenvalues of a time-dependent potential 
at a fixed time. They therefore shift with time. However, in order to study the effects 
of this shifting in more detail one needs to solve the Bethe-Salpeter equationliHil which 
describes the dynamics of the correlated electron-hole pairs. This is beyond the scope of 
our model calculation which is intended to proof and explain the existence of transient 
excitonic states in Cu. In the T — K direction exists a sp-like band which crosses the Fermi 
energy. In addition to these single-particle states there are also a wealth of excitonic states. 
Assuming that the lifetime of excitonic states is longer than the electronic lifetime of band 
states (see Sec. [V]) this explains the experimental finding that the lifetime of electrons in 
these two directions are very similar. It also explains why GW calculations - which only 
address electron-electron scattering of excited band states with valence states - reveal shorter 
lifetimes as compared with experiment. For both directions the eigenvalues were calculated 
using the total potential t>tot( r ; t) at t — 0.01 fs. In the T — X direction the excitonic energies 
are a little bit higher due to the high energies of the lowest unoccupied band (7.5 eV at the 
X point). The excitonic energies in this direction were calculated at t — 0.002 fs. Due to 
the large slope of the lowest unoccupied bands the transient excitonic states which can be 
probed by experiment are all in the vicinity of the high symmetry points X, L, and K. 
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V. LIFETIME OF ELECTRONS IN TRANSIENT EXCITONIC STATES 



Having shown the possibility of transient excitons in metals we want to close with some 
comments about the lifetime of electrons in these states. The dominating process which 
determines the lifetime of excited electrons in band states is electron-electron scattering. 
Due to collisions with electrons of the valence band the excited electron is scattered out of 
its state, creating an electron-hole pair. This process - which diagram is shown in Fig. [3[a) 
- and exchange scattering are the two processes which are considered in GW calculations. 
However, only diagram ^.a) contributes to a final lifetime. Experience shows that the lifetime 
due to this process is proportional to the density of states of the final states. In metals this 
is therefore a rather efficient process. The bound states as described by transient excitons 
on the other hand can only decay via pair annihilation - pair creation, i.e., via an Auger 
process. The decay process is displayed in Fig. ||.b). As a consequence of the bound-state 
character of electrons trapped in a transient excitonic state we conclude that the lifetime of 
these states is longer than in band states. This reasoning would explain the disagreement 
between ab initio GW calculations for the lifetime of excited electrons along the V — K 
direction and the experimental data. 



VI. CONCLUSIONS 

In this article we presented a model which allowed us to calculate transient excitonic 
states in metals. The model was applied to jellium and to real Cu. The results obtained 
within the jellium model describe how a (nearly) free electron system screens a localized 
positive charge. An analysis of the time-dependent total potential shows that the complete 
screening of the external charge takes 10 to 20 fs. This is the time needed until the polar- 
izability is fully built up and in equilibrium. At this time the potential is nearly identical 
to the Thomas-Fermi potential. This near identity is caused by the fact that the screening 
process is completely determined by the plasmon excitations of the system. The contribu- 
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tions from single-particle transitions - which are always present in the response function for 
finite wave vectors - are negligible. The potential possesses bound states for the first 0.2 fs. 

The situation in Cu is different since the imaginary part of the density-response func- 
tion (which is proportional to the loss function) has no distinct plasmon peak. Instead 
the density-response function of Cu shows plasmon-like excitations originating from the d 
electrons. This double peak sits on top of a broad structure caused by single-particle tran- 
sitions. The screening of a localized charge distribution is therefore faster. Nevertheless, 
the total potential possesses bound states for small times. In our model calculation for the 
determination of excitonic energies these eigenvalues are plotted with respect to the lowest 
unoccupied band. We find that the excitonic energies are in the energetic region of states 
detected in the measurements of the lifetime of hot electrons in crystalline Cu using TR- 
2PPE spectroscopy. It is important to stress that only experiments with very short laser 
pulses allow the detection of these states. 

Our calculations explain thus the puzzling experimental fact of states which are no band 
states. We interpret them as transient excitons. With the coarse estimate of the lifetime of 
these states we are also able to explain the difference between the experimentally determined 
lifetime of hot electrons in the (110) direction with respect to the results of GW calculations; 
we think that the states probed in this direction are not band states but excitonic states. 

Our model does not yet include the possibility to determine the exact lifetime of electrons 
in transient excitonic states. The model is only applicable up to the moment at which the 
electron is trapped in a bound state and an exciton is formed. Work along this line, i.e., the 
evaluation of diagram b) in Fig. ^] is in progress. 
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FIGURES 
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FIG. 1. The band structure of Cu calculated using a plane-wave code£j For the energy range 
between 1 and 4 eV above Ep Cu possesses no band states in the V — L and the V — X directions. 
There is a sp-like band in the V — K direction. 
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FIG. 2. The first four non-zero components of the hole potential v h (r) calculated from the 
density distribution of a d electron at the X point in Cu. The solid line is the spherical sym- 
metric component (1=0, m=0). The dotted, dashed, and dot-dashed lines denote the components 
(1=2, m=0), (1=2, m=l), and (1=4, m=4), respectively. The long-dashed line shows the Coulomb 
potential multiplied by \/4vr. 




q/k F 

FIG. 3. The total potential v to t(q,t) for r s = 2.07 in atomic units. The upper panel shows 
vtot(q, i) for q = 0.1 &f plotted as a function of time for different values of rj in a small time window. 
The broadenings used are 0.01 eV, 0.1 eV, and 0.5 eV (solid, dotted, and dashed line, respectively). 
The oscillations with the frequency u p can be clearly seen. In the lower panel Vtot(q,t) is plotted 
as a function of q in units of the Fermi wave vector for t = 10 fs. 
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FIG. 4. Snapshots of the total potential vt t(r,t) as a function of r for different times. The 
upper row shows the potential for small times calculated using r s = 5. The horizontal lines denote 
the eigenvalues of the potential calculated with a reduced mass of ix = 0.5. In the lower row 
vtot(f>£) is shown for larger times and r s = 2.07. It can be nicely seen how the system overscreens 
the perturbing potential which results in a decaying oscillation around the asymptotic form of the 
total potential. Note that the times shown cannot be used as a measure for the frequency of the 
oscillations which is much larger (see Fig. ||). The dashed lines denote the Thomas- Fermi potentials 
for the two Wigner-Seitz radii. 
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FIG. 5. Influence of the vertex correction / xc on the total potential. The upper panel shows 
Vtot(q, t) calculated with and without the vertex correction (dot-dashed and solid line, respectively). 
In the lower panel v to t(r,t) is displayed. The solid line shows again the result of the calculation 
with / xc = and the dot-dashed line the result of a calculation with a vertex correction. In both 
plots the dashed line denotes the Thomas-Fermi potential. The calculations where done for r s = 5 
and t = 1000 fs. 
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FIG. 6. The imaginary part of the density-response function xg=G'( c Li w ) of Cu calculated for 
three wave vectors q + G. The underlying calculation of the polarizability was done including 200 
bands. 
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FIG. 7. Four snapshots of the total potential in Cu. The total potential was calculated using 
a polarizability which includes 200 bands. Note the different scales in the upper and lower row. 
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FIG. 8. Band structure of Cu in the three main crystallographic directions (solid lines). In 
addition to the single-particle states the transient excitonic states at t = 0.01 fs (L — V and V — K 
directions) and t = 0.002 fs (T — X direction) are plotted as horizontal bars. 
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FIG. 9. Diagrams for the possible decay mechanism of electrons in a band state and a transient 
excitonic state. In a) the diagram for the decay of a band state is shown. This diagram determines 
the lifetime of excited electrons in GW calculations. Diagram b) is the annihilation process involved 
in the decay of a transient excitonic state. 



24 



